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TKCUNK'AL MKMOUANDUM 7HMH 


STATUS REPORT: DATA MANAGEMENT PROGRAM 
ALGORITHM EVALUATION ACTIVITY A1 
MARSHALL SPACE FLIGHT CENTER 

INTRODUCTION 


It in widely recognized tiuit the moat prolific data neneratora have lH*en 
multis|>eetral Ima^in^ senaors, aueh aa theae uaed on I^ndaat, and the Synthetic 
A|H*rture Itadar, aueh aa thoae to Ik‘ uaed «>n Si>aaat. Theae typea of aonaors 
generate a large volume of data at a high rate, which ofUm requlrea complex 
analyala procetluix's to convert the data to uaahle Information. Thua, tlio 
magnitude and complexity of the data ofUm cix>aU> costly bottlenecks In ground 
data processing syater<.s. To overcome this pix)hlem, atandanl cost-effectlvo 
preprocessing and processing prt)cedurea that pi*oduce accurate and consistent 
results are needed; tht-se do not exist at present. The difficulties In developing 
these procedures have I'»ecn determined to be causi'd by the natuix’ and form of 
the data, as well ns by tlu‘ approaches ludiig taken to its analysis. 

Ihuler the Data Management I’logram, an algorlUim evaluation activity 
was Initiated to study the processing prol)lem by assessing the indeperuk'nt and 
Interdependent effects of registration, compression, and classification U'ch- 
niques on landsat data for several discipline applications. The objective of the 
activity was to make recommendations on selected applicable Image |)rwesslng 
algorithms in terms of accuracy, cost, and timeliness or to propose other 
solutions to the processing problem. As a means of accomplishing this objec- 
tive, an Image Ctxllng Panel was established with members from several NASA 
Centers who represented different ty|)es of Image processing Interests. The 
Panel established the evaluation criteria, selected thi* data sets and algorithms, 
performed the analyses, and presented the results during the Panel meeting. 
AfU'r each meeting, a panel summai'y' was written which reported on the |)cr- 
formanco of Individual algorithms, combinations of algorithms, or on new 
recommendations that were made. 


APPROACH/ PLANNED ACTIVITIES 


Under approved HTOP 050-21-02-02, there la a requirement to (l) 
acquire and maintain detailed inforuTUtion on all technique development in the 
agency and elsewiv,'rc, with emphasis on OA s|K)nsorcfl activites, esj>eclally 
within the Data Management and Karth Resources Programs; (2) analyze, 
evaluate, and l)rinn to focus information on technique development; and (3) 
report and advise Headquarters on emphasis chanj>es and the need for technolo;©' 
development. In partial fulfillment of these requli*ements, the f)bjectlves of 
this effort are to update and document all source lma};e processing programs, 
to determine the comhined and i nde inuidtMit effects of registration and com- 
pression on classification, and to evaluate n'j>istration, compression, classi- 
fication, and clian;;c detection !*esults. 

I)»;rin}; UY77, the fliscipline applications of land use inventoi^’ and 
aurlculture were to ixj studied usin^, six I-andsat passes over Mol)ilo Hay, 
Alabama, for land use and eleven I^indsat passes over a IJ\CIE ( lar^e Area 
Crop Inventory Experiment) supersite ( Einncy County, Kansas) for agriculture. 
The particular processing techniques to l)c studied In conjunction with these test 
sites w'erc as follows: 

a. Renistratlf)n 'rechniques: 

1. Nearest Neiulilior (NN) 

2. Bilinear Interpolation (HI.) 

3. Hlcul)lc Literpolatlon (HC). 

p'or a brief description of the rcfristration techniques, see NASA TM X-73348, 
"Nearest Neighbor, Bilinear Interpolation and Bicubic InteriTolatlon Gconraphlc 
Correction Effects on I.;mdsat Imaninjj, " R. Jayroe, Marshall Space EllKht 
Center, Septeml»er 1975. 

I). Compression 'Pechnlques: 

1. Adaptive Differential IHilse CckIc Modulation (ADPCM) 

2. Adaptive llaearmand (All) 

3. Combination AH/ ADPCM. 
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Tor a brief description of tliese techniques, see TilW Final Keport No. 20.®3G0, 
”Studj" of Adaptive Methods for Data Compression on Multispectral Scanner 
Data,” Contract NAS2-SG94, NASA Ames Research Center, March 1977. 


4, Cluster Coding Algorithm (CCA) — See ”A Joint Clustering/ Data 
Compression Concept,” Ph.D. Diesertation, University of Sou thei’n Califbrnla, 
E. E, Hilbert of Pasadena, Jet Propulsion Laboratory, May 1975. 


H. G. Mulk, 


Blob Algorithm — See ”Uscrs Guide for the Blob Program Package,” 
Code 9bG, Goddard Space Flight Center, June 23, 1977. 


e. Classification Techniques: 

1. Gaussian Maximum Likelihood (kLlXLK) 

2. Sequential Linear Classifier (SLIC) . 


For a brief description of these techniques, see NASA TN B-S240, "Classification 
Software Technique Assessjuent,” R. Jaj-'roe, et al. , Marshall Space Flight 
Center, Mayl97G. 

3. Gaussian Mixture Density Classifier (GMDC) 

4. Iterative Clustering Classifier (ICC). 


For a description of these programs see "Users Guide of Software Documenta- 
tion for the Algorithm Simulation Test and Evaluation Program,” J. Lyons, 
Johnson Space Center. 


The preceding image processing techniques were chosen because they 
are representative of all the known approaches used in I’egistering, compressing, 
and classifying multispectral imagery, and they appear to be generally accept- 
able to analysts. The specific tasks using these techniques to be accomplished 
bv MSFC within the Image Coding Panel in FY77 are described in the foil 
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paragraphs. These tasks are being worked by General Electric conti’aetor per- 
sonnel under contract NASS-324u4, as ’well as MSEC in-house personneL 


The emphasis of the majority of the tasks is to determine what has 
changed in the data, how much change has occurred, and where has the change 
taken place as a result of registration and/or eompression. Secondly, how' are 
these diang'es manifested in the classification results The dominant factor in 


the evaluations is the claMsirieution accurncy foHf)wo(l by the eost effectiveness 
of the analysis procedure^ with the* result l>cln^» that there is usually a tradeoff 
bet^^•eon the two. 


Task 1: Updating, Modification, and Documentation 
of Software Techniques Library 

A 8oft^^•are techniques library had l)ocn established at MSFC under the 
Data Management Proj;ram prior to the formation of the Image C<xllng Panel. 
This library is augmented with image processing techniques that are utilized by 
the Panel and facilitates the transfer of these techniques to the I’anel members. 
The source programs for the software are housed on a disk pack under the name 
SMTOOl, are In I'ortran, and are run on an IP^I-.^^»()/ 7r». The task Is to docu- 
ment the programs In *'reatlablo Fortran,” prepare a software techniques manual 
of all programs and a user’s manual for updating and utilizing the library. 
”lloadable Fortran” is a (k)cumcntatlon concept that requires maximum usage of 
comment statements that arc Indented in a U)glcal fashion to explain the software 
coding and such that the program can Ik* read In a manner similar to a textbcK)k. 


Task?: Determination of Registration 
Effects on Image Data 

Perform NN, PL, PC' geographic correction on the IJVt'IK and Mobile 
Pay test sites and save the corrected data sets. To (k'termlnc and document the 
registration effects, use t!ie following approach: 

1. Compare channel histograms, means, and variances for the original 
and geographically corrected Image test sites 

2. Compare joint channel histogi*am pairs for the three registration 
techniques and the original data 

3. For Identical channels, compute image difference histograms and an 
actual Image showing where the differences occur for the three registration 
techniques. 
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Task 3; Determination of Compression 
Effects on Image Data 

l^oiir eomprossioii teeliniqutis will be yeleetocl from a eombinaHoii of 
those previously meiitloiietl. For the selected toctmlques and three compression 
ratios, compress, reconstruct, and save the image data for the two test sites 
and use the following ap|)roach: 

1. Conipare chaniicd histograms, means, and variances for the original 
and compressed/ reconstructed image test sites 

2. Compare joint channel histogram jjalrs for the four compression 
techniques and the original data 

d. For identical diannels, compute image difference histograms and an 
actual image showdng where differences occur for the compressed/ reconstructed 
and original images. 


Task 4: Determination of the Combined Effects 
of Compression and Registration on Image Data 

Select two data sets from each test site tliat have been geographically 
corrected using NN, BL, and EC. Select the corresponding two data sets that 
have been compressed at three different ratios and reconstructed using the four 
Compression techniques. GGograijhieally correct the compressed/ reconstructed 
images using NN, BL, and BC, and compjress and reconstruct the previously 
geographically corrected images using the four comin*ession tecimiques and 
three comijression ratios. The resulting imag;os will be saved for future use 
and differ only in the order of compression and registration. To determine 
and document the combined effects of image correction and compression, the 
following approach should be used foi’ the two gxmips of images; 

1. Compare cliannel histograms, ivseans, and variances for both gx’oups 
of images and with the original images 

2. Compare joint channel histogram pairs ter both groups of images and 
wdth the original images 
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8, For identical eliacnela, comimte image 
actual image allowing wJiero dl£ft?reneea oeetir fort 


3rence hiatogTams and an 


a. The two graupa of image a 

b. The geographically corrected — compressed/ reeonstrueted 
images and the geographically corrected images of tasiv 2 

c. The coiiipresaed/ reconstructed — geographically corrected 
[mages and the compressed/ recoiistructea images of task ii* 

Task 5; Determination of the Effects of Compression 
and Registration on Classification 

Using the elassificatitm programs Mi\XLlK, SLIC, MI>:ture Density, 
Iterative Clustering Classifier, and the digital ground truth maps (GT'^l) for 
both test sites, classify the image data and comi»are pr-';el by pfKel with the 
ground trudi data to obtain tlie following classification accuracies for evaluation; 

1. Classification of each original image set compared with tlie GTM 
geographically distorted to overtaj- v;ith original data 

2, Classification of each geograpiiically corrected (NK, BT, BC) image 
set compared wdth the GTM 

a. Classification of each compressed/ reconstructed (four lecluniques, 

3 compression ratios) image set compared with the GTM geographically distorted 
to overlay with the eompressed/ reconstructed data 

4. Classification of the image sets described in task 4 compared with 
the GTM 


5. Glassification/ GTM difference maps of the above to show wlie^’C 
disagTeemen ts occur. 
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Task 6; Determination of Change Detection Effects 

From tht* I^CIK and Mobile Hay imaKea, aelcct the KOOKraphieally 
corrected data acta (ualntt NN, HI., or HC) that ap|>eur beat auited for dctectlnj; 
acaaonal or other tvtx>a of chanKca baaed upon the roaulta of the previoua taaka. 
For each teat alte and correapondlnn aeaaonal ImaKoa, nao the followlan a|)|)roach 
to determine change detection cffecta: 

1. To determine htm’ the data chan^ea with aeaann, compute difference 
hlato^rama arid joint hlatonrama, and Ima^^c diflerencea l)otwoen correapondInK 
lmaf;e banda 

2. Compute claaalflcatlon map difforencea and contingency tablea to 
determine how' the claaalflcatlon reaulta change with acaaonal data 

3. Uae the CITM aa a maak to extract feature Information from the 
aeaaonal data (including and ixjludltiK tranaltlon pixela) by computing joint 
hlatn(;rama between corrc'apondlnn banda to determine how thia Information 
mlnratca through the data apace aa a function of acaaon 

4. Uae the claaalflcatlon mapa aa a maak to extract Information from 
the (lata (Including and excluding tranaltlon pixela) by computing joint hlato(;rama 
beHveen correapondlnj; banda and c(jmparlnn with the ground truth feature 
extracted Information 

T). For the compreaaed and Keonraphleally corrected data, |>erform the 
comparlaon 1 to 3 of thIa taak to determine the cffecta of compreaalon on change 
detection. 

Aa a rcault of the Imapie CcHlInn Panel activities, four additional rcKjul ce- 
ments have lx»on identified to assist In the transfer and evaluation of linage data 
and ima(;c analysis products. These requirements consist mainly of the develop- 
ment of software programs to measure linage entropy, to utili/.e fractional pixel 
accuracy In uround truth/ classification map assessments, to facilitate Imaso 
data and analysis product transfer, and to explore aiid evaluate a compression/ 
classification procedure In terms of analysis eost reduction potential. The 
tasks to Ixj accomplished for these requirements are discussed In the followlnR 
para^^rapha. 
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Task 7: Measurement of Imaye Entropy 

Devolo|> and document In **roatlabIe l orlron" a aoftwaro program to tie 
included in the MSFC Sottwarc l ilirary tliut measurcH the entnifiy of an ImaKC 
aa a function of lipnd. The foIlowlnK is a veriiul description of the desired soft* 
v.are proKt'am. 

Differential entropy is often used as a measure of ima((o information or 
activity. It actually ix'presents the minimum data rate tvquired for encodinK 
differences of the ituBKe with a /ero-memory enctnier, under the assumption 
that the differences liave a stationary distriinition timHi^hout tiie imaKC. In 
Keneral, the differences an« the ik'viations iietweon an estimate for each picture 
element (liased on all pn’ccdlnn picture elements) and tlic picture element itself. 
In practice, it is adequate to use as an estimate for eacli picture element the 
averajjc of the element aliove and ix'hind; l.o., let 




A 
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be the esMmate for X 


Ij’ 


where* 
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X 




Note: Where I xj Is the nearest inU'jter value to x (ties may either 
way). Then the ellfferencc for element (i, j) is 
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Clearly, by aoiullnn a MtartltiK relerenco picture element, the rest of the Imayc 
can be tranjmltted exactly In terms of the dlffercncea. Di'flno as the fre- 
quency of cK'currenc** of a difference having value ft averaged over the entire 
ima^o. Then the two-dlmet>slonal difference entropy for that Imane Is 


M . - X I'j lOK, . 

ft 

The ontrt>|)y II represents the lowest bit rate possible for encoding the differ- 
ences. Thus II cun Im' a ust'ful measurt' of Um' difficulty one would h£.ve In 
ccxllnj; a particular imuKc, assuming an exact reconstruction of the ImaRc is 
desired. Kntropy measures are to be calculated for <110 Mobile Hay and 1J\CIE 
data sets. 


Task 8: Data Transfer Format 

t)evelop and document In ”ri ailablo Fortran’* a sj^ftwaro program to be 
Inc’- dr i n the MSFC Software l ibrary for linage data transfer that has the 
^luow'lnn output characteristics: 

1. Data tape should Im> 7 track, r>5(t bpi atui contain data within the ran^e 
/.ero (dark) to 2b." (bright). 

2. Size of Lest site should be approximately 2.'>() pixels by 2.*)r) records. 

3. There should be one band of ima);e data |)or file and onj record of 
(lata per scan, with each flic having a header record. The header record 
should contain loo characters or loss with K bits per character. 

•1. Training areas are to l>e polygons with at most nine vertices, and 
Interior |x>lnts are to l)e utilized as tralnlii}; data. If a ti*alnlnj; area has only 
tw'o vertices, a imarest neighbor rule is used for solcctinn training data. 

T). A listing of the training area vertices and class identification should 
accompany the data tape along with a histogram of each band of data. 

<5. For pixel location, each image starts at pixel one and reconI one, 
rather than (o. 0 ) . 
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1 % A liytiog of tlie gi^ogi’ajjhie eoiTcction ti’aiiafui*matl6« should also bo 
inelttded with tlw data tajju. 

Imago data tapes and analysis products are to be made available for transfer 
reciuest. 


Task 9t Fractional Pixel Accuracy 

Opdato and doeumeiU. In ’’readable Fortran” a software irrogram for 
inclusion In the MSFC H»:s.tlwaro library to utiHae Iraetlunal piKel accuracy in 
the ground truth/elassineatiun ma[j {CITM/CM) aescssments. The approach 
utilised in the soliwarc developmc!ii corresponds to digitising the GTM at a 
higher fesoluilon than tlie CM such that gt*veral neurcst neighbor pjlKels on the 
GTM can be used to evaluate one piriel on the CM. The updated software evalua* 
tion prograiii has to allow for different CTM resolutions and for assignment of 
statistical weights used with the GTM nearost neighbor idnels. This approach 
is to be evaluated using tlie procedure described in NAflA TM “73347, ’’Evalua- 
tion Criteria for Bodwaro Classification Inventories, Accuracies and Maps,” 
by R. dayroe and compared with the roistilts obtained from the classification 
a. 8 ses 3 mcnt using the Mobile Bay and LAFIE data sets. 


Task 10: Reduced Vector Representation 


This task is a research effort direcced toward reducing the cost of 
analysing multispc"‘C’al imagory. When four band, six bit image data are 
acquired, the number of distirict or different vectors possible is (04)“’ =* 

16 777 210. For Landsat data, this maximum possible Is limited to 7 081 6UU 
different feature vectors, which is the number of picture elements contained 
In one Landsat frame. The question arises, ”Ilow many different feature 
vectors are required to represent an itnage and still satisfy all user require- 
ments V” If tlie immber of different feature vectors cuiild be reduced to 204S, 
for example, the number of classification eomr>utations could be reduced by a 
factor of Several thr'tusand and result in considcjrable analysis cost savings. 

The foHovvdng aijpi’oach is recommended to evaluate the potential cost savings; 

i. Develop a sotlware program for computing tlie number and population 
of different feature vectors contained in an image. 
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7. A listing of the KC«Kra|jhic correction transforniati(Mi 8hr>uld al8o Imj 
included with the data tape. 

linage data tapes and analyais pixxlucts are to lx> made avallal)le for transfer 
upcjn request. 

Task 9: Fractional Pixel Accuracy 

Tpilhtc and document in ”rea<lal)le I'ortran** a software pro>;i*am for 
inclusion in the MS7C Software I.il»»ary t<» iitili/o fractional pixel accurac*' in 
the KnAiiul ti*uth/ classification map (CJTM/t'M) assessments. 'I’he appix>ach 
utilized in the software developn>ent corresp<*nds t<» diKiti/in^; the CITM at a 
hinher resolution than tlu* CM such that several nc*!irest neiuhlmr pixels on the 
GTM can be used to evaluate c*-.e pixel on the C'M. The updated s«jftware evalua- 
tion proKram has to allow for different fi I’M resolutions and for asslunment of 
statistical weiuhts used with the (ITM nearest neinhbor pixels. This approach 
is to 1 h' evaluated usinn the proceduiv described in NASA 'I'M -7iJ;i-17, *’Kvalua- 
tion CriU^ria for Software Classification Inventories, Accuracies and Maps," 
by H. Jayroe and compared with the results obtained from the classification 
assessment usin^; the Mobile bay and lA'TFJ data sets. 

Task 10: Reduced Vector Representation 

I'his task is a research effort directed toward rt*duclnf> the cost of 
analyzing multlspe< i ’al imagery. When four band, six Ijit lmano data are 
acquired, the number of distinct or different vectors possible is (bl)^ = 

1() 777 21t). For Landsat data, this maximum possible is limited to 7 r)Hl 000 
different feature vectors, which is the number of picture elements contained 
In one Landsat frame, 'rhe question arises, "How many tliffercnt feature 
vectors are required to represent an Imu^fc and still satisfy all user require- 
ments?" If the number of different featui’c vectors coidd be reduced to 2048, 
for example, the number of classification computations could bo reduced by a 
factor of several th tusand and result in considerable analysis cost savinjis. 

The fc>llt>win>;' api)i*oach is recommended to evaluate the potential cost savings: 

1. IX'velop a softwai’c prof»ram for com|)ulinn the number and po|;ulation 
cf different feature vectors contained in an imaj;e. 
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2. Utilize a version of Ukj cIusU*r coding ainorithm (CCA) for redue'ng 
each M by N pixel array to a maximum of 1(» different feature vectors, and 
develop a program for meruinn tin* resulting numiier of feature vectors con- 
tained in the imaKc down to a desired final numiMir (i.e. , the final number is 
an input variaido to the Pi'onram). This will l)o an operational program that 
provides an output tai)e containing a tal)le of the final dosirt'd numlwr of differ- 
ent fei tui*c vectors, their po|)ulations and lalx'Is, followed i)y a table of their 
ima^e locations that arc identified liy tlieir lal>els. 

3. IX'velop a software program for convertinK the output tajxs back to 
feature vector format. 

i. Develop a software classiricatlon program that utili/cs the above 
output tape (step 2) lor automatic classification in the /following manner: 

a. A user selects training areas and desi^^nates one class numl)cr 
for all feature vectors contained in the training area, ilms, the classification 
program will have a table that relates a class lal)el to the feature vector laixjls 
and an itnage Iwation. By usinn tlx* approach in C’CA, additional feature vectors 
bclonninn to a particular class can also be Idetnifit'd and labeled. 

b. The user will have an option of producing; a map and population 
talde contalninn only those feature vectors identified in the training; procechire 
as belonniuK to a particular class, or of invoking a classification procedure 
that assigns all of the feature vectors to one of tlie selected classes. The 
unassigned or unclassified feature vectors can l)c classified in the following 
manner. Compute the centroid of tiu* fealuix' vectors belonging to each class 
and store be result. Tor an unclassified feature vector, find the closest two 
centroids and then compare the unclassified feature vecU)r witii individual 
memlxji'S or representatives of the two identified classes. Assign the feature 
vector to the class of the closest individual member. The user should also 
have tlie option of outputting a one cliannel CM or a classification feature vector 
map using the centroids. The classT.'cation results ol)tained from this method 
should be evaluated on the 1J\CIK and Mol)ile Bay data sets and compared with 
the other techniques using the |)ix)ceduix> previously descriljed. llie effects of 
linearly stretching the data before applying the feature vector rcch:ction proce- 
dure should also be Investigateti together with the cost effectiveness of the 
procedure. 
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I)nrinK 1<T7K, tlie applicutionu of mineral and |>etroleuin rcsoureea and 
water and environmental quality will Im* Htudied in a similar manner. The 
al^o..ithm evaluation aetivity is ex|K»eted to Ik* eoncluded l)y FY7U, 


ACCOMPLISHMENTS 


Durinn tlio period Oetf>l)er 15)7(1 to September 1977, two Imai;c Coding 
I*anel meetings wore ludd and the teehnical results obtained l)y tl>c panel mem- 
bers were documentt'd in llic November 30, 1970, and IIk* March 23, 1977, 

Image CcMling Panel Summarit*s. 'I'lie summaries deserlbed and discussed the 
various independent and Interdependent effects »)f registration, compression, 
and classification on tlie iJald Knol), Tennessee, quadrangle teat site. The 
most significant technical accomplishment was a recommendation and j)ro|K)sed 
prwedure lor using a reduced vector representation to approximate Landsnt 
multi S|)ectral Imagery. This recommendation provided the capability of reducing 
l^ndsat processing costs by tlm*e orders of magnitude and of compressing the 
data l)y a factor of tme to thrc*e orders of magnitude. This recommendation was 
further exploix^d by developing sf)ftware to determine the number of different 
vectors in an imago and, as a result, to determine a measun^ of image com- 
plexity. Image complexity was then used as a criterion for evaluating registra- 
tion and compression techniques to determine their capal)llity with the reduced 
vector representation. 'I'he results of this study were pi’osented a'. NASA 
Headquarters on Septeml)er 7, 1977, and are currently in |Hil)licatien. 

At tile request of Headquarters, an assessment of tin* Image Coding 
Panel’ s accomplishments was fonvaivled to Headquarters bv the Panel chairman 
on Juno 29, 1977. The- assessment addressed wliat ha. I lieen accomplished, the 
usefulness of the work, the results that were obtained, prolilem areas and 
possilile solutions. 

Specific accomplishments relative to tlie previously descrilied tasks under 
planned activities are listed in fhe following paragraplis by month These 
accomplishments addressed the activities that are iieing conducted at MSFC, 
but the Panel participation and assistance in analysis and evaluation provided 
by the Jet Propulsion lalioratory (JPI.) and the Ames, (loddard (GSPC), 

Johnson (JSC), and I^ingley (lallC) Space Cenb'rs should be recognized and 
acknowledged. Although this rec(;f,nition can be liest aclileved liy reviewing the 
Panel summaries, some of the cooiieratlve efforts can be evidenced in the 
progress that Is reported. 
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May 1977 

1. ProtfroHa . 

a. liocation of Ground Control Points (GCP) (Task 2). Tho OlKital 
linage Itectlfiration %stom (I)IUS) was requested and received from GSFC. 

The program has the (•apal)lllty of ^cneratin); printer plots of GCP nolfjhhorhoods, 
which are manulflcd by cubic lnter|)olatlon, and of removlnn the effects of Earth 
rotation and sensor sampling delay. It was found that this process removed 
significant distortions In the region of the first Mobile Hay GCP. These dis- 
tortions had piwlously caused difficulty in correlating tlx* Imagery with quad- 
rangle map^. However, the DIRS prt)gram for this stop stores the output Image 
In core and hence requli*ed 27-1 K bytes of storage. Consequently, a program to 
perform these corrections by lines was written and tesled using l^andsat Imagery 
of tho Huntsville Jetport main runway. The Improvement in ap|)oarancc was 
striking. DIRS was ix>quested to Improve MSEC’ s registration ca|>ability In 
prt)vidlng geographically corrected test site data to tho Image Coding Panel for 
subsequent analysis and documentation. 

b. Registiatlon ('lask 2). During tho process of rectification, 

DIRS removes systematic effects such as oversampling, along line Earth rota- 
tion, sensor delay, mirr<»r velocity nonlinearity, and Earth curvature. How- 
ever, it df>es not partition the image to rotate segments of the Imago. Conse- 
quently, tho imago is not rotated, and the core requirements arc large (580 K 
bytes). The program does allow specification of a rotation angle, which could 
bo accommcHlated by jiartRlonlng the imago. 

c. Data Compression ('I'ask 3). The programs supplied by NASA 
Ames via TRW (the adaptive compression techniques) liave been analyzed and 
appear to be readily amcnalde to modification to use the s|wclficd data sets. 
Additional supporting programs received from TRW i)y mall were successfully 
compiled. 

d. Additional Software Analysis (Task l). Tho DIRS package Included 
specialize Fortran callable input-output routines used on GSFC* s niM-300. These 
routines are being compared with wliat is normally used on MSI-'C* s DiM-360 to 
determine what would be more efficient at MSFC. 
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0 , Mlseella-ueoufcs. At GSFC’ b recjuost the d^ekB iov the HINDU 
elaaeiOctitton pi^ogram (Iluerlatleally Inai/u'ed Neighborhood Discerning' 
Unsiiperswed) and a source tape oi' TIIW’ a cum|')ression j-^rograms v/ere checked, 
updated, verified, and loaded on tape tor shipment to GSFC. 


2, Proijlenis. The DIRS program abnormally ends during transformation 
of the image corners to Universal Transverse Mercator (UTM) coordinates. A 
deck has been setup to insert debug siatements into specified subroutines in 
hupos of Isolating and klentifying the problem. 


June 1977 


1. Progress. 


a. Updating, Modification, anti Docmnentation of Bofhvare Teclmique.s 
Library (Task 1 ). The registration software package received from GSFC and 
the data compression software received from TRW have been analysed with 
respect to flo jlity to the stated techiikiuea and relationships used In the algo- 
rithms. The registration sofUrare, as received, would require dedicated use 
of the 360 computer due to core requirenieiits. The compression software 
requires some g'eneraliS4atlon to be used with data sets of vaiy/lng siaes. A new 
library of routines which can be used to satisfy externa! references during link- 
age editing has been created, and is named IMAGE. ANALYSIS. LIBRARY. 


b. Determination of Hegistration Effects on Image Data (Task 2) . 
The problem v/hlch arose In modifying certain subroutines was due to incon- 
sistent subroutine names in the source programs supplied, and was resolved 
satisfactorily. The two programs which rectify images with and without rotation 
were successfully executed on a small test area and on a Mobile Bay Test Site, 
■respectively. The IBM double precision least squares routine, DLLSO, was 
used to obtain registration mapping polynomials of degrees one through five. 
Some aspects of stating the sizes of the fit errors obtained and effects of fit 
aeeiumcy on input data magnitudes (due to computation errors) are under 
consideration. 


e. Detei’niination of Compi’cssiori Effects on Image Data (Task 3). 
The final briefing of the TRW study on data compression was held by telephone 
conference on June 22 and was attended by the GE personnel assigned to this 
contract. The NASA centers represented were Headquarters, Ames, JSC, 
MSEC, and GSFC. 
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o. MI«c‘clIan<‘oua. At (ISFC* a r<*quoHt tho dockM for the IIINUU 
claHHiflcatiun program ( Ilucri»tie:illy Inspired Nci^hlxirhotMi DiscorninK 
rnsujx*rawed) aiul u .source tu|x* of TUW* h compression programs were checked, 
updated, verified, and loaded on ta|)c for shipment to (ISl C. 

2. l*rol)lems. The DIKS pn>Kram al)normally ends during; transformation 
of the imane corners to I'niver.sal Transverse Mercator (TTM) cwirdinates. A 
deck has been set up to inse^'t debu^ statements into s|X‘cified subroutines in 
ho|x*s of isolating and identifying tlK> problem. 


June 1977 


1. Prf>u ress . 

a. UiKlatin^, Modification, and Uocu.nentation of Software Techniques 
Library (Task 1). The registration software package received from CiSFC and 
the data compression software received from THW have lx»en analyzed with 
respect to fli...llty to the stated techniques and tvlationships u.setl In the al^o- 
rithms. The registration software, as received, would require dedicated use 

of the :i(i() computer due to core requirements. The compression software 
requires sf)me neneralization to be used with data sets «>f varying sizes. A new 
library of routines which can be u.sed to satisfy external references during link- 
age editing has been created, and is named IMAGK. ANALYSIS. LilfHAUY. 

b. DeU^rmlnation of Ite^istration Effects on Ima>j;e Data (Task 2). 

The problem which arose in modifylnf* certain subroutines was due to incon- 
sistent subroutine names In the source programs supi>lied, and was resolved 
satisfactorily. The two |)roM;rams which rectify images with and without rotation 
were successfully executed on a small U'st area and on a .Mobile liay Test Site, 
respectively. The IBM double precision least .squares routine, DLI-.SO, was 
used to obtain registration map|)ln{; polynomials of decrees one throuj;h five. 
Some aspects of stating the sizes of the fit errors obtained and effects of fit 
accuracy on input data mannitudcs (due to computation errors) are under 
consideration. 


c. Determination c>f Compression Effects on Imane Data (Task 3). 
The final briefinn of the TUW study on data compression was held by telephone 
conference on .June 22 and was atUMuled by the (IE |)crsonncl assinned to this 
contract. The NASA centers repre.sented were lloaflquarters, Ames, JSC, 
MSEC, and CSFC. 
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2. Probleina. There ap|)ear to bt* son.'* IneonsUtenclos In the techniques 
for correction of Karth rotation, sensor sampling, and mirror motion In the 
reulstratlon program. There aiv two versions of «)ne of tin* data compression 
source programs. AH(' should verify whic-h version Is correct. 


July 1977 

1. IM’on ress . 

a. r|Klatln) 4 , Modification, and l>»cumentation of Software Techniques 
Mhrary (Task l). ('onversion to r(>adahle Kortran has been accomplished for 
tht* three data compression techniques recommemk’d In the TItW study. These 
are block adaptive transform (Cosine or iiadamard), bbu-k Adaptive Differential 
Ihdse Code M«)dulatlon (ADl’CM), ami a hybrid technique. Most of the routines 
have iK’on placed in the linage Analysis l.ibrary in loa<l module form. 

b. Determination of Iteulstration Effects on Imane Data (Task 2). 

The problem with least squares filtinn to obtain the rc^;istration mapping poly- 
nomials was due to the larue values of the CTM Input eoordinates. The routine 
appears to operate satisfaetorlly if the means or the minima are subtracted from 
the coordinates. 

The Inconsistencies In the DlltS received from C.SKC were acknowledfiod 
by Its author to be errors. 

A niethod for partitioning U) apply DIRS to lar^e imaj;es has Ijcen for- 
mulated. It la planned to generate output I mane senments based on the nrid 
cells defined by tlx* <‘oarse mappinn functi«)n nrid. To facilitate this process, 
the beninninn and endinn output pixels and eastinn coordinates for each cell have 
been placed In tables, rather than beinn ivcomputed. 

c. Determination of Compression Effeets on Imanc Data (Task 3). 

The updated data compression software is beinn validated by repetition of the 
orlnlnal checkout runs. Preliminary results have Ix'en obtained from the applica- 
tion of the ADPCM alnorithm to I.ACIE data. The compression operation ap|)cars 
to quadruple the numlwr of unique feature vectors In a 1J\CIE sample senment. 
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(1. MiacoDaiu'outi. Nonadaptivo data c'ompn'tiHlon alKorithrna, which 
were delivered to MSFC i)y I'llW, are lK*lnj{ pivpared i»y MSFC for tranamittal 
to (ISFC. Tlie adaptive U*cht)iqiies wen' sent previously. 

The di(>lti/.ed (ITM of Mol)ile Bay is bcinn examined at MSFC in an attempt 
to find sources of ern)r, such as highways which are «ilfflcu1t to di((iti/.e and are 
u(’casslona1ly missed. 

2. Prol)lems. Hequlred moiiiflcations and corrections to the reRistratlon 
prop;ram a»*e Krt'ater than anticipaU'd for handlitiK a test siU* as larjte as Mobile 
Bay. 


August 1977 


1 . Brourt'ss . 

a. l'|M^atin^^, Mcxilflcation, and Documentation of Software Techniques 
l ibrary ( Task 1). Documentation (*f three liata compression teciaiiques in 
readable Fortran is undei’way. 'i’he techniques are i)lock adaptive transform 
usin^ a Cosine (TRANC) or Iladamard (TRANH) transform, i>lock ADPC'M, and 
a hvl)rid comi)inatlon of the al>ove two techniques ( if^’BltDC-ADPCM/ COSINE 
and IIYBRDH-ADPCM/ lladarmard) . 

b. Determination of lU'^istration Effects on linage Data (Task 2). 

The OSFC DIRS has been mo<lified so that lar^e imanes ran l)e correcU'ci in 
segments. The segments are inised on tiu* >;rid cells defined by tlie coarse 
mapping function grid. I’nder tlx' constraint of the memory availaldo for holding 
input Itnagc ri'cords, the numl)er of output map cells that can i)c generated in 
each row is computed. If any row cannot Ih> completed, the image is divirk'd into 
vertical segments, geometrically corr(*cU'<l, and the output line segments are 
then merged. 

c. IX'termination of Compression Effects <*n Imago Data (Task 3). 
iXiring validation of the daUi compression software, impro|)cr handling of 
floating point to integer round-off was discovered in the Hybrid and ADPCM 
programs. Correction I’csulU'd in an approximately 25 jn'roent deervase in the 
Mean Squared Error ( MSE) between the original and reconstructed images. 


I 


I 
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I 

i 
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The three data coiniiresaion techniquea have lx»en applied to a 1 j\CIL’ 
sample seKnient. and the numlM>r of unique feature vectors in the reconstructed 
imbues uas computi'd. This is a measure of the linage activity which has been 
ad<led l)y the compression/ reconstruction process. The averane bit rate required 
to transmit the compressed ima^^es was specified as 2-bits/ pixel. The MSE 
and numlx»rs of feature vectois are as follows: 


Method 

MSE 

No. Veci 

2-1) Iladamard 

2.09 

5925 

ADI'CM 

1.85 

7200 

Hybrid-Had. 

2.07 

7r):i.3 

Original 


0552 


Ta|H?8 of the original and reccnstructed images were ol)talned for use on 
the Test lied film record(*r. 

d. I)(‘terminatlon of the Effects of Compression and He^istration on 
Classification (Task 4). The ciustcrinjj program liLOIi was received in a 
virtual memory version which is not avaiialde at MSEC. The projjram has been 
modified and is executing on MSEC*s IIiM-.'5<iO/ 7r>. (’ontour and cluster maps 
for 1*ACIE data have i)een ol)tained. 

e. Eractionai I’lxei Accuracy In Ground Truth/ Classification Map 
Assessment (Task 9). Avalialiie documentation on diKitai manipulation of GTM 
and assignment of class niunlK'rs is l)einn studied as l)acknround for deveiopinK 
a soft\vare proj;ram to use the GTM di^iti/ed at a hipher resolution than CM. 
Since the software program will use the procedure described in NASA 

TM X-7.‘i.'M7, ’’Evaluation Criteria for Software Classification Inventories, 
Accuracies and Maps,” an exercise In familiarization with this procedure is In 
progress. 

f. Entropy Measure (Task 7). The Image entropy measure soft- 
ware has been programmed and debugged. 

% 

g. Reduced Image Vector Representation (Task lO). Software for 
computing the number and population of different feature vectors contained In an 
image has been prtigi’ammed and (iebugged. 
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h. MidccUancmis. UVCIt and Bald Knob Ijmdaat teat aitc data 
have been tranamilted to Jl*l. and JSC, reapoctlvcly. JSC ia providing MSKC 
with the Ciauaaian mixture denaity and Iterative cluater elaaalflera and ia 
runnini; a teat eaae on tht> Bald Kn<»b data. 

2. Pr«»blcM))8. The converalon and ch'.eklnn of the program KI.I.TAB 
waa jMjatponed In anticipation of reeelvlnn an updaU'd jr»0 veraion. 

Extensive uae of the Teat Bed waa not planned becauae of disk hardware 
prol)1cma. 


September 1977 


1 . Proureaa . 

a. Updatin{{, Modification, and Documentation of Software Techniques 
I.ll)rary. The Bl.t)It pro^jjram package pn>vlded l)y CSKC waa modified to run on 
a machine without virtual memory. Tlie routines have l)een placed in MSEC* a 
Imanc Analysis I.ll)rary in load mo<hde form. Tin* adaptive data compression 
techniques (transform, differential pulse c<kIc nuKlulation, and a hybrid tech- 
nique) Imve been documenUMi in readal)le Fortran, 

1). IX^termlnation of Bet^istration Effects on linage Data. The dis- 
tortions in l>tindsat MSS imagery caused by sensor delay anti Earth rotation are 
quite pronounced in mannlfi(‘d imagery and should Ik* removed prior to selecting 
pround control points. A propram was written to extract tiM; data sumuindlnp 
a candidate pround control point location, renun'e the ciistortlons, and mapnify 
usinp cubic interpolation, 

c. Determination of Compressitm Effects ott Imape Data. D:ita 
compression was applied to a lJ\(TK sample scpmc*'t and the December 1973 
pass coverlnp Mol)lle Bay. The ivsults were analyzed l»y computinp moan values, 
histoprama, difference imapes, MSE, and the numix'r ol unique feature vectors 
in the Imapes. Numerical results are piven In Tal)lo 1. 

d. Fractional Pixel Accuracy in Gmund Tnith/ Classification Ma| 
Assessment. To dipiti/e the 1.AC1E (ITM, a 15 l)y 13 in. trtinsparency of the 
map (supplied i>y JSC) was reduced to a 2 liy 1.8 in. nepative l)y the MSFC 

p hotop r a ph 1 c 1 abo r a to ry . 


Method 


MSK 


S|K'cd 

( Pixels/ sec) 


No. 

Vectors 


_1 ^ L 

Ij\CIK Sample Se>;ment 


Original 



3 552 

ADPCM 1 Bit 

7. 92 

Hi.'} 

10 492 

ADPCM 2 Bits 

1 . sr> 

050 

7 206 

2-D Badamard 1 Bit 

4.89 

457 

5 985 

2-1) Badamard 2 Bits 

2.09 

299 

5 925 

2-1) Cosine 1 Bit 

9.71 

147 

5 290 

2-1) Cosine 2 Bits 

1 . 02 


5 754 

Bybrld-Badamard 1 Bit 

7.91 

59(i 

G 896 

Bybrid-Badamard 2 Bits 

2.07 

504 

7 033 

Ilybrld-Coslnc 1 Bit 

f).01 

280 

7 101 

Hybrid-Cosine 2 Bits 

1.94 


7 684 

BLOB® 

21.90 

209 

558 


Mobile Hay Data 


OrlKlnal 



27 090 

ADPCM 1 Bit 

4.87 

040 

>58 000 

ADPCM 2 Bits 

1.14 

024 

40 845 

2-1) Badamard 1 Bit 

2.47 

459 

37 870 

2-1) Badamard 2 Bits 

1 . 19 

957 

35 406 

Bybrld-Badamard 1 Bit 

5. 01 

740 

41 248 

Bybrid-Badamard 2 Bits 

1 . 90 

547 

46 575 











This uegutivo was digitized by a mierodenslometer (seamier) at se: 
Intervals of fie, 2d, and 12. S g. These seaiinlng intervals reimesent resolatlons 
(magnifications) ot 6X, lOK, and 2uX, respectively, relative to Landsat data, 
Tlio oatput of the scanner ^vaa stored on three tapes, one for each of tlie throe 
seatmlng intervals. 


To facilitate plotting of the I^ACIE data, the 50 g scanner tape was copied 
to another tape in Fortran format using an IBIVI yjfstem/hdo general I/O packap. 


A Line iTiiiter plot was made of tlie map using tlie Fortran tape and the 
IBM 8ho eomfjtiter, This plot will be used lO obtain tlie coordinates of ground 
control points for I’oglstratlon of the GT.M^ witli the CM. 

As an aid in soleethig; training samples for elaysifying data uii the Landsat 
tape, the data on the GTM were ealsaified manually; i.e. , tabulated data of the 
GTM wex’e transferred to a photocopy of the map. 


All existing cempiiter program (SLlC) for linear el a; 
Landsat data was modified to accommodate the LACIF data. 


3ifieatlon of the 
A computer pro» 


gram supplied by GSFC for geometric, registration of the GTM waa reviev/ed. 
This program will be used to fit (determine piKol correspondence of) the GTM 
to CM. 


October 1977 


1. Progress . 

a. Updating, Modification, and Uocumetitation of Software Techniques 
Library. The histogram routine being used In the analysis of compression 
results was modified to show the occurrences of negative values in the difference 
images. 


The algorithm (SLIC) used for dctermiumg linear classification functions 
was modified to store intermediate results, resulting in a computation time 
saving. The maximum likelihood classification algorithm was modified to 
include thresholding. 


This iiCKailvc was dlKlti/od by a ndcnKlonslomoter (scannor) at scanninf; 
intorvals of TiO, 2o, and 12.5^. These seannInK Intervals ifprewnt resolutions 
(maKnIfirations) ol f>X, lOX, and 20X, resiwetlvely, ridatlve to iiindsat data. 

TIh? output of the seanner was stored tm thive tapes, one b»r each of tht? three 
seannlnt; Intervals. 

To facilitate plotting of the IJVCIK »lata, the r>0 p scanner tape was copied 
to another ta|K? In Fortran format using an IBM system/ .’U>0 general I/O package. 

A I. Inc Printer plot was made of the map using the Fortran tape and the 
IBM atiO computer. This plot will be used lo (»blaln tiu* c<K)rdlnaU's of ground 
control |X)lnts for registration of the GTM with the CM. 

Ah an aid In selecting training samples for classifying data on the landsat 
tajK!, the data on the (ITM were calssified manually; I.e. , tabulated d?ta of the 
GTM were transferred to a phoBwopy of the map. 

An existing computer program (SI.l(’) for linear classification of the 
l^ndsat data was imxllfled to uccomnuxlate the data. A computer pro- 

gram supplied by CiSFC for geometric registration of the (ITM was reviewed. 

This program will Ik' used to fit (determine pixel correspondence of) the Cl’l'M 
to CM. 


October 1977 


1 . Progress . 

a. Updating, Modification, and Documentation of Software Techniques 
I.il)rary. The histogram routine Ix'ing used in the analysis of c impression 
results was modified to show the occurrences of negative values in the difference 
images. 

The algorithm (SI IC) used for <ktermhong linear classification functions 
was mcxllflcd to store Intermediate results, resulting in a computation time 
saving. The maximum likelihcxxl classifir ation algorithm was modified to 
include thresholding. 
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b. Dt'tcrmination of Rf>((i8tration Effects on ImaKe Data. The 
ground control |K>ints for tlk* lk*ct*mlj<‘r 73 pass of Mobile Hay data were pre- 
viously selecte<i usintt a niat;nifieation of -IX. The Kn^^ind control |M>int coor- 
dinates were rodeterminofi usin^ a magnification of 10, and with removal of 
distortions caused l>y sensor delay and Earth rotation. lx*ast square fits were 
determined for |)olynomlals of deurees one through five. The HMS fit ern*rs 
for the corrected Imagery wen* lower by an average of 2H |K*rcont. The HMS 
errors are nlven in Table 


TAHI.E 2. HMS EHHOHS IN ITHTNC (IHOl'NI) CONTHOE POINTS 


Polynomial 

Denree 

With 

Distortion 

Without 

Distortion 

1 

0.H97 

0.G19 

2 

0. H57 

0. r>7H 

3 

0.7ii7 

(». 5.37 

•1 

0.()45 

0.4G0 

5 

0.319 

0.405 


c. Determination of Compression Effects on linage Data. Additional 
data compression results are ulven in Table 3. llisto>;rams of tlie difference 
images (orininal — reconstructed) for ‘I baiuls (A 1.AC1K data uslnn ADPCM at 

3 bits/pixel are niven in Flj^ire 1. IMots of MSE and number of feature vectors 
In the reconstnicted imaf^es are nlven in Eij^^ires 2 and 3. 

d. Fractional Pixel Accuracy in Clnnind Tnith/ Classification Map 
Assessment. landsat data of the IJXCIE site were received as files of a 
maprictic tape with each file correspondinj; to a particular day that the data were 
obtained. Several files of this tape wore copied to a disk to simplify accessing 
the data. 

A Line Printer Plot was made of File 9 (Scanner Data of May fi, 197fi), 
and training; samples, 20(» each of wheat, summer fallow, corn, sornhum, and 
^rass, were selected carefully to lie within the field boundaries of the plot. 
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TABLE 3. DATA COMI’KESSION EFFECTS 




No. 

S|>eed 

Total 

Mf»h<Ml 

MSE 

Vectors 

( Pixels/ sec) 

CPC Time 

1J\C1E Sample Senment 


OrlKlnal 




■1 

ADPCM 1-Dlt 

7.92 


0.30 

HBI 

2-Dits 

I . H.3 


0.30 

HiB 

3-Dits 

0. 710 


007 

mSM 

2-1) iladamard 1-Hit 

•1.H9 

.3 38.3 


0:47 

2-Hits 

2.09 

3 923 


1:03 

.3-Hits 

0.930 

3 893 


1:14 

2-1) Cosine 1-Hlt 

3.71 

3 230 

147 


2-Hits 

1 . 02 

3 7.34 

141 


.3-Hits 

0.0.39 

.3 800 

124 

2:. 33 

Hvbrld-!la<lamard 1-Hlt 

7.31 

0 830 

.390 

0:30 

2-Hlts 

2.07 

7 033 

304 

0:43 

3-Hits 

0. 77.3 

0 299 

.310 

0:42 

Hybrid-Cosine 1-Hit 

.3.01 

7 101 

280 

0:77 

2-Hits 

1.91 

7 081 

238 

1:31 

3-Hlts 

0.71H 

0 308 

222 

1 :37 

1 hi.oh‘^ 

21.3 

.3.38 

209 

1:00 

Mobile Hay Data 


1 

Original 
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ADPCM 1-Hit 

■1.87 

78 .334 

r,40 

37:10 

2-Hlts 

l.M 

40. 84.3 

024 

38:28 

2-1) Iladamard 1-Hlt 

2. 17 

37 87<1 

4.39 

.32:19 

2-Hlts 

1.19 

33 400 

337 

07:08 

Ilybrid-lladamard 1-Hlt 

3.01 

41 248 

740 

32:20 

2-Hits 

1.90 

4(J .373 

.317 

43:31 

Fast Karhunon-Loeve 

•1.11 

48 2.33 

100 

145:00 

i-nit*> 






a. Applied to area of 112 x 112 pixels as opposed to 1 12 x 192. 

b. Not |<jrformed unck'r this contract. 
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OF POOH QUALITY 
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P'lRuro !• Difforenco iinaRCS (original — reconstructed) for 4 bands 
of data usinj; AI)I’('M at .'{-bits/ pixel. 


MEAN SQUARE ERROR 





File 1 (Scannor Data of ()ctolK.‘r 22, 1975) was clussifiefi usin^ a maxi- 
mum likelihood algorithm t»> sr *arato wheat Inim the oth<;r crops. It was found 
that wheat represented 40 iwreent of the crops Krowln^ at that time. 

The HAMTFK-1*1)1*-1 1/ 15 Imane Processlnn Facility was used to remove 
annotations from the dlultl/ed GTM of tlu* l-AC'IE data, and a polarold print was 
made of the map to verify that all of the annotations liad been removed. 

The etiltcd GTM was stored on a tape In Fortran format for Keometric 
correction and reulstratlon at a later <late. 


RECOMMENDATIONS AND CONCLUSIONS 


The Panel should be continued throunh FY78 and possibly Uirounh part of 
FY79 to satisfactorily compiete the work that has been initiated. The Panel has 
now attained a position where the evaluations can l)c conducted in a routine 
manner and contributions can be made toward determining; the best and most cost 
effective pi*ocedures for reulsterlnt;. compressing, and classifylnj; I^ndsat 
lmaf;cr>'. The time consuminn and seemingly unprcxiuctive efforts of ncttliif; 
orjtanizod, develo|)lnj; an analysis plan, establishing; evaluation criteria, etc., 
were one time activities that have Ix'cn substantially completed and, hence, 
more effort can be coi.ccntraU*d rui the analysis and evaluation. 

Two milestones have been achieved by the I'anel which indicate that sub- 
stantial |)i‘oy[rcss is beint; made. First, there is sufficient technical maU?rial 
boint; Kont^rated and presented by the participatinj; cenUn’s to justify extending; 
the originally scheduled one day I’anel nu'ctiiiKS to two days. Secondly, an 
issue has surfaced In the technical results that wei*e recently presented, which 
place the l^anel at a crossroad with rcf;ards to pursuing certain types of rej;istra- 
tion and compi'ossion techniques. This issue is centered ai’ound the number of 
unique four dimensional vectors contained in a I^mdsat multispectral imase, 
which can also be used to achieve a rt'duced processing; cost by uslnj; a hlstof;ram 
typo format. 

The histoRram format Is obtained by cxtractinp; all of the unique vectors 
and tlToir number of occurrences fr«)in a liindsat imaRe, or a portion there of, 
placInR that information at the boRinninR of a data tape, and followinR it with a 
description of the iinaRc. The linaRc is described at each picture clement 
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i 

I 
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location i)y one number which iiiontifica the unique four dimensional vector that 
l>eIonKS thew^. Important cost retiuctions for most ty|ics of processing can be 
realized, by procesaitif; each unique vector once and '^pplyinR the results to the 
picture clement locations witii a tai)le look-up procedure, since many of the 
unique vectors can occur imndreds or thousands of tintcs in an ima^e. 

An issue arises with the rcRistration and compression techniques that 
utilize spatial avcraninn or an equivalent, ix?causc they tend to generate large 
numbers of new unique vectors from the original data i)y interpolation or 
requantization. This issue l)ocomes more important if significant cost reductions 
are to Ik* achieved by appi'oximating an image with a reduced numl)er of unique 
vectors. For this case, i-oglslration and compression tec*uiiques must have the 
property of not generating any new unique vectors to preserve the processing 
cost reductions. 

Thei'e arc t >\’0 main |K)lnts in using a hlstogi*am format that need to be 
emphasized. First, the histogram format does provide a certain amount of 
information preserving compression, depending on the imago data (typically 
1.5:1). Secondly, the format directly illustrates tliat what is l)eing processed 
and interpreted in a multlspoctral image is in actuality the unique vectors and 
not cverj' picture element, which represent a processing overhead when unique 
vectors are repeattHi in an image. In this context, the function of a picture 
clement is only to show wlierc tin* unique vectors art' located in the image scene. 
Thus, instead of i>ing prt)cesslng costs to the numl)er of picture elements in an 
Image as is commonly done, rethiced costs can l)c achieved l)y using a histogram 
format and tying the costs directly to the elements that are actually being 
analyzed, i.c. , the number of unique vectors. One aix?a that needs to l)e explored 
is the use of information -preserving eompressien techniques (which typically 
achieve compression ratif>s of two or tlirec to one) on the liistogram formatted 
data. With the combination of techniques, it may Ik? possll)le to achieve informa- 
tion preserving compression ratios of three or five to one. 

If the approach of lying the processing costs to the numlK?r of unique 
vectors is extcndci. to n'duce costs even further, then the numl)cr of unique 
vectors in an image has to lx? reduced. This reqidres that an imago be approxi- 
mated l>y a reduced vector representation, which in reality is a classification 
map with a relatively largo number of classes. A final classification map with 
8 or 10 classes, for example, is a reduced vector rcpi’osentation carried to an 
extreme. Two ajx?as of inquiry that need to lie pursued arc the tradeoff between 
reducing the prtjcessing costs by reducing the numl)er of unique vectors versus 
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allowing an analyst Uie flexibility of conducting whatever type investigation is 
desired, and determining the ix*st proceciurc for ol)taining a reduced vector 
representation tliat can Ixj automatically implemented. One approach has l)cen 
develo|XHi, and is currently under investigation, for providing a reduced vector 
representation. The approach is an extended version of the Cluster Coding 
Algorithm (CCA) called the C'ascaded CCA. The CCA combines groups of 
similar unique vectors into cluster centers, wldle the Cascaded CC'A combines 
groups of similar cluster centc'rs producing a significantly reduced vector 
representation. Iloth of these software programs arc compression schemes, 
which allow the user to process the data in compressed form. The concept of 
using a reduced representation appears reasonable when it is recognized that a 
largo image contains lens of thousands of unique vcctf)rs and to 50 |)crc’ent 
of these vectors only occur once. If the tlcslred final result is a CM containing 
less than 100 classes, for example, then the large numl)or of unique vectors 
appear to represent an over abundance of variation of information. Hy approxi- 
mately reducing the number of unique vectors in a Landsat image to 2000 in 
number, it ap|x?ars |>ossibic to achieve processing cost rcihiction factors of 
approximately .*1000 and noninformation preserving compression ratios of 10:100. 
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STATUS REPORT: DATA MANAGEMENT PROGRAM 
ALGORITHM EVALUATION ACTIVITY AT 
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By Holwil H. Ja;, roc, Jr. 

'I’hc Information in this report lias Ik'cmi rcvif»wcfl for security classi- 
fication. Itcvlow of any information concerning Department ol Defense or 
Atomic Enerny Commissica proj^rams has lieen made l»y tla* MSIT Security 
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unclassified. 

This document has also lK*en reviewed and approved for technical 
accuracy. 



